Reply to "Comment on 'Evidence for the immobile bipolaron formation in the 
paramagnetic state of the magnetoresistive manganites"' 
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Zhao et al. (Phys. Rev. B 62, Rll 949 (2000)) reported studies of the oxygen-isotope effects on 
the intrinsic resistivity and thermoelectric power in the paramagnetic state of several ferromagnetic 
manganites. The isotope effects on the intrinsic electrical transport properties are not consistent 
with a simple small-polaron hopping mechanism, but can be well explained by a model where 
there coexist immobile bipolarons and thermally-excited small polarons. Recently, Banerjee et al. 
(Phys. Rev. B 68, 186401 (2003)) wrote a Comment on our paper, where they misunderstood our 
theoretical model, erroneously thinking that we explain the electrical transport in the paramagnetic 
state as due to bipolaron hopping. They confuse large polarons with small bipolarons, use incorrect 
conditions for bipolaron formation, and even incorrectly reproduce our data. We show that the 
model used in the Comment is inconsistent with any features of the observed oxygen-isotope effects. 
More evidence is provided to support the theoretical model used in our original paper. 



Experimental evidence for small polaron charge carri- 
ers in the paramagnetic state was provided by transport 
measurements [1]. It was found that the activation en- 
ergy Ep deduced from the conductivity data is one order 
of magnitude larger than the activation energy ob- 
tained from the thermoelectric power data. Such a large 
difference in the activation energies is the hallmark of the 
small-polaron hopping conduction [1]. In our recent pa- 
per [2], we reported studies of the oxygen-isotope effects 
on the intrinsic resistivity and thermoelectric power in 
the paramagnetic state of several ferromagnetic mangan- 
ites. The isotope effects on the intrinsic electrical trans- 
port properties can be well explained by a model where 
there coexist immobile bipolarons and thermally-excited 
small polarons. In our theoretical model, the electrical 
transport in the paramagnetic state is due to small po- 
laron hoping, but such polarons are thermally excited 
from localized bipolaronic states. Recently, Banerjee et 
al. [3] wrote a Comment on our paper, where they misun- 
derstood our theoretical model, erroneously thinking that 
we explain the electrical transport in the paramagnetic 
state as due to bipolaron hopping. They confuse large 
polarons with small bipolarons, use incorrect conditions 
for bipolaron formation, and even incorrectly reproduce 
our data. We show that the model used in the Comment 
is inconsistent with any features of the observed oxygen- 
isotope effects. More evidence is provided to support the 
theoretical model used in our original paper [2]. 

Fig. 1 shows the temperature dependence of the 
resistivity of the oxygen-isotope exchanged films of 
Lao.75Cao.25Mn03 and Ndo.vSro.sMnOa. The data are 
the same as those reported in our original paper [2] . The 
only difference is that the resistivity is plotted as a func- 
tion of T /Tc in our original paper [2], where Tc is the 
Curie temperature. Comparing the data shown here and 
the data shown in Fig. 1 of the Comment [3], we can 
clearly see that they do not reproduce our data correctly. 
Further, the authors of the Comment [3] fit their incor- 



rectly reproduced data by 

p = DTeMEp/kBT), 



(1) 



where the coefficient D is inversely proportional to the 
characteristic phonon frequency z^o- 
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FIG. 1. The temperature dependence of the resistivity of 
the oxygen-isotope exchanged films of a) Lao.75Cao.25Mn03 
and b) Ndo.ySro.sMnOs. The data are the same as those 
reported in Ref. [2]. 



1 



We can also fit the data of Fig. 1 by Eq. 1. The curves 
fitted to the resistivity data of the oxygen-isotope ex- 
cliangcd Lao.75Cao.25Mn03 films are shown in Fig. 2. 
The best fits yield the fitting parameter D = 7.36xl0~^ 
ricm/K for the ^^^O sample and D = e.lSxlO""^ ficm/K 
for the ^^O sample. We can see that the D value for the 
^^O sample is larger than that for the ^*^0 sample by 20%. 
Since D oc l/i'o oc a/M, we expect that the D value for 
the ^^O sample should be smaller than that for the ^^O 
sample. This is in sharp contrast with the isotope effect 
on the D value obtained above. This implies that Eq. 1 
is not a relevant formula to describe the electrical trans- 
port in the paramagnetic state of Lao.75Cao.25Mn03. If 
the authors of the Comment would provide the fitting pa- 
rameters for both isotope samples, they would not claim 
that their model could better explain our data. 
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original paper [2], we fit the resistivity data of the ^^O 
and samples of Lao.75Cao.25Mn03 by Eq. 2 and find 
that the parameter C for the ^^O sample is larger than 
that for the ^*^0 sample by 34%, in agreement with the 
theoretical model (Eq. 2). In fact, this model can quan- 
titatively explain [2] the combined isotope effects on C 
and on Es- 

On the other hand, if small polarons are bound to lo- 
calized impurity states rather than to localized bipolaron 
states, Es would be independent of the isotope mass in 
the zero-order approximation according to Austin and 
Mott [6]. In our original paper, we have only consid- 
ered this possibility [2]. Here, we shall consider the other 
possibility, that is, a polaron is bound to an impurity 
site with the screened Coulombic interaction reduced by 
a factor of the static dielectric constant (continuous 
medium approximation) . The binding energy of the im- 
purity state is given by [7] 



13.6 m* 

Ec = eV, 



(3) 



where m* is the effective mass of a polaron and rUe is the 

mass of a free electron. The thermally excited polaron 
density from the impurity states is [7] 

n = {2Na)y^{m*kBT/2T,h^f/^ exp{-Ej2kBT), (4) 

where Nd is the dopant concentration. The thermoelec- 
tric power is correspondingly given by [6] 



S='^{EJkBT + a'), 

e 



(5) 



FIG. 2. The temperature dependence of the resistivity 
of the oxygen-isotope excliangcd films of Lao.75Ca().25Mn03. 
The solid linos arc the curves best fitted by Eq. 1. Tlic fitting 
parameter D for the ^^O sample is larger than that for the 
sample by 20%. 

In fact, we have shown [2] that if fesT is less than the 
polaron bandwidth Wp and only thermally excited po- 
larons from localized bipolaron states contribute to the 
electrical transport in the paramagnetic state, the resis- 
tivity is given by 



p=-^exp{Ep/kBT), 



(2) 



where Ep = Ea + Es, C = (aVe^v^)(l-05Wp)i-V/ii'o, 
a is the hopping distance, Es is the activation energy for 
polarons to be thermally excited from localized bipolaron 
states, and Ea depends on the polaron binding energy 



where Es = E^jl and o! is a constant depending on 
the kinetic energy of the polarons and on the polaron 
density [6]. With the mobility of polarons [8] /U oc 
Vo/kBT exp{—Ea/kBT), we finally have 



p = BT"''eMEp/kBT), 

^3/4 



(6) 



where Ep = E^ + Ec/2 and B oc Wp' /vo- 

In Fig. 3 we fit the resistivity data of the isotope ex- 
changed Lao.75Cao.25Mn03 by Eq. 6. The best fits yield 
B = 1.13x10"^ ncm/KO-25 for the ^''O sample and B 
= 0.89x10"'' f7cm/K°-2^ for the ^^O sample. It is clear 
that the B value for the sample is larger than that 
for the ^^O sample by 27%. This is in agreement with the 
theoretical prediction. If we assume that is inversely 
proportional to the square root of the reduced mass of the 
manganese and oxygen atoms, Vo will decrease by 4.6% 
upon replacing by ^^O. Using the oxygen-isotope ef- 

3/4 ^ ■ 



Ep and bare hoping integral t. The quantity C should ^ct on B and the relation B Wp we find that 

Wp for the *'0 sample is smaller than that for the °0 



strongly depend on the isotope mass M and decrease with 
increasing M. This is because Wp decreases strongly with 
increasing M according to Wp = 12t cxp(—TEp/huJo) = 
12texp(— g^), where F is a constant (<1) [4,5]. In our 



sample by 46%. The exponent of the oxygen-isotope ef- 
fect on Wp is Po = -d\nWp/d\nMo = 3.2, where Mq 
is the oxygen mass. 
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FIG. 3. The temperature dependence of the resistivity 
of the oxygen-isotope exchanged films of Lao.75Cao.25Mn03. 
The soUd Unes are the curves best fitted by Eq. 6. The fitting 
parameter B for the sample is larger than that for the 
1*0 sample by 27%. 

From the thcrniooloctric power data of the isotope ex- 
changed Lao.75Cao.25Mn03, we have found that [2] Eg = 
13.2 meV for the ^""O sample and Eg = 18.7 meV for the 
sample. Then we have E^ = 26.4 meV for the ^'^O 
sample and Ec = 37 A meV for the ^^O sample. Using 
Eq. 3 and Wp oc 1/m*, we find Po ~ 3.0, close to the 
value (3.2) deduced from the resistivity data above. This 
appears to indicate that this model can quantitatively 
explain the isotope effects. 

From Wp = 12texp{—rEp/?iu)o) = 12fexp(— 5^), we 
can easily show that (3o = g^/2.6. Here we have used the 
fact that LOo is inversely proportional to the square root 
of the reduced mass of the manganese and oxygen atoms. 
Using i3o = 3.0-3.2, one has g"^ = 7.8-8.3. Assuming that 
the bare effective mass of charged carriers is about 0.5me, 
we find that the effective mass of polarons m* is about 
1500me, which seems unlikely. 

It is worthy of noting that Eq. 3 is valid only if the Bohr 
radius rs = 0.53esme/m* A is significantly larger than 
the interatom distance (2 A) [7]. Using = 20 (Ref. [11]) 
and m* = ISOOrrie, we get = 0.007 A and Ec = 51 cV. 
The calculated ii^c = 51 eV is over three orders of mag- 
nitude larger than the measured one. The calculated tb 
= 0.007 A implies that the continuous medium approx- 
imation does not hold. Thus Eq. 3 cannot be applied 
to doped manganites unless the effective mass of charged 
carriers is close to me- This is possible only if the charged 
carriers are not small polarons, in contradiction with the 
large isotope exponent /3o and the observed small po- 
laron hopping conduction. Therefore, the model where 
polarons are bound to impurity sites with the screened 
Coulombic interaction (continuous medium approxima- 
tion) cannot explain the observed isotope effects on the 
intrinsic electrical transport properties. 

On the other hand, if polarons are bound into much 



heavier bipolarons [9], the random potentials produced 
by dopants can completely localize bipolaron band. The 
minimum separation between the polaron and bipolaron 
bands is A = 2(1 - r)Ep - Vc - Wp, where Vc is the 
Coulombic repulsion between bound polarons [5]. It is 
clear that the value of A could be less than 0.1 eV, in 
agreement with the measured value of about 30 meV in 
Lao.75Cao.25Mn03. 

Within the same model where polarons are bound into 
the localized bipolaron states [2], we have also deter- 
mined the bandwidths for the isotope exchanged samples 
of Lao.75Cao.25Mn03, that is, Wp = 49.0 meV for ^^O 
sample and Wp = 38.8 meV for the '^^O sample. The val- 
ues of Wp lead to /3o = 2.0 and g"^ = 4.6. Then the bare 
hopping integral t is calculated to be 0.41 eV from the 
values of and Wp. The inferred value of t is in reason- 
able agreement with the band structure calculation [12]. 
Further, since E^ = A/2 and A = 2(1 - T)Ep -Vc- Wp, 
the change of Eg due to the oxygen isotope substitution 
is SEg = 6A/2 = -6Wp/2. With dWp/2 = - 11.2 meV 
for Lao.75Cao.25Mn03, one has SEg = 5.5 meV, in quanti- 
tative agreement with the value (5.5±0.6 meV) deduced 
independently from the thermoelectric power data [2] . 

The electron-phonon coupling constant g'^ can be also 
determined from optical conductivity data. It is shown 
[5] that the optical conductivity contributed from pola- 
ronic charged carriers exhibits a broad peak feature at an 
energy of Em — 2g'^TuOo. Jimg et al. [10] have identified 
the polaronic conductivity and determined the value of 
Em as a function of doping x in Lai_a;Caj,Mn03 system. 
From their results [10], one can easily see that Em ~ 0.7 
eV for X = 0.25. If we take fiojo = 0.074 eV (Ref. [2]), we 
have g'^ = 4.7, in quantitative agreement with the value 
(4.6) deduced from the isotope effects. 

In our data analyses [2] , we assume that bipolarons are 
completely localized by random potential due to dopants 
and defects, and the conductivity is only contributed 
from thermally excited small polarons. The authors of 
the Comment misunderstand our theoretical model, erro- 
neously thinking that we explain the conductivity in the 
paramagnetic state as due to bipolaron hopping. Fur- 
thermore, they confuse large polarons with small bipo- 
larons. From the value of a' in Eq. 5, which is found to be 
less than 1 for our samples, the authors of the Comment 
conclude that [3] our data are consistent with small po- 
laron hopping conduction mechanism, but disagree with 
bipolaron hopping conduction mechanism. This is in- 
correct. The value of a' can only discriminate between 
small polaron and large polaron conduction mechanisms. 
Bipolarons are not equivalent to large polarons. 

Finally, the arguments [3] against the formation of 
bipolarons in doped manganites do not have scientific 
ground. The authors of the Comment argue that bipo- 
larons cannot be formed in doped manganites because 
there is not enough disorder. We do not believe that 
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disorder is a required condition for the formation of 
bipolarons. Bipolarons can be formed when electron- 
phonon interactions are strong enongh to overcome direct 
Coiilombic interaction between two polarons [5] . The au- 
thors of the Comment [3] also mistakenly think that the 
charge disproportionality is a necessary condition for the 
formation of bipolarons. Intersite bipolarons do not in- 
volve any charge disproportionality. In this Comment, 
there is also a mistake in the criterion for the adiabatic or 
nonadiabatic hoping conduction mechanisms. The bare 
bandwidth should be used to compare the parameter (f) = 
{kBTEp/TTf/'^{nuJo/'!Tf''^. But the authors of the Com- 
ment [3] use the effective polaron bandwidth instead of 
the bare bandwidth to compare the parameter </>. It is 
clear that the bare bandwidth is larger than (f) so that the 
adiabatic hopping conduction is relevant in doped man- 
ganites. 
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